Chirped quasi-phase-matching with Gauss sums for production of biphotons 



(N 

o 



(N 

S 1 : 

Oh- 
^— > ! 

c : 

=5 

cr 



> 
On 

in 

m 

in 
o 



X 



D. A. Antonosyan, 1 ' 2 ^ A. R. Tamazyan, 1 ' 2 '0 and G. Yu. Kryuchkyan 1 ' 2 '!! 

1 Yerevan State University, Alex Manoogian 1, 0025, Yerevan, Armenia 
2 Institute for Physical Researches, National Academy of Sciences, 
Ashtarak-2, 0203, Ashtarak, Armenia 

We study the theory of linearly chirped biphoton wave-packets produced in two basic quasi- 
phase -matching configurations: chirped photonic-like crystals and aperiodically poled crystals. The 
novelty is that these structures are considered as definite assembles of nonlinear layers that leads 
to detailed description of spontaneous parametric down-conversion (SPDC) processes through the 
discrete Gauss sums. We demonstrate that biphoton spectra for chirped photonic crystals involving 
a small number of layers consist from definite well-resolved spectral lines. We also discuss the 
forming of broadband spectra of signal (idler) waves in SPDC for both configurations as number of 
layers increases as well as in dependence of chirping parameters . 

PACS numbers: 42.65.Lm, 42.50.Dv, 42.65.Yj 



I. INTRODUCTION 

One of the challenges in the fields of quantum op- 
tics and quantum communications is generation of new 
sources of nonclassical light, such as biphotons with con- 
trollable spectral and temporal properties. The standard 
method for the generation of biphotons is spontaneous 
parametric down-conversion (SPDC) [3-[|, where photon 
pairs are generated under action of a strong pump field 
interacting with a nonlinear crystal. In particular, of con- 
siderable interest is preparation of biphoton wave-packets 
with small correlation times between two-photons as well 
as with broad frequency spectra. Correlated photon pairs 
with broad spectra are mainly generated by SPDC using 
method of quasi-phase-matching (QPM) that provides a 
feasible alternative to conventional phase matching for 
many optical parametric process applications. Several 
methods have been suggested and experimentally imple- 
mented in this direction for generation of biphotons, in- 
cluding SPDC in aperiodically poled crystals [J-|9(. Stud- 
ies in this direction allow us to understand three-photon 
interactions at scale of the shortest possible correlation 
time that is a single optical cycle [8]. Thus, these inves- 
tigations are interesting from the fundamental point of 
view as well as are useful for further applications in quan- 
tum technologies. QPM leading to photon pair produc- 
tion is effectively realized in layered nonlinear structures, 
particularly, in periodically-polled second-order crystals 
or photonic crystals. The applications have been done for 
synthesis of twin photon states by manipulating overall 
group delay mismatches between interacting waves in a 
multilayered structures by using compensating dispersive 
effects |10l - |l3| . for production of entangled three-photon 
states in cascaded parametric processes [T3 - [To1 |. 

As it has been recently experimentally demonstrated 
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biphotons with broad spectra and ultrashort correlation 
time at a high emission rate can be generated using 
chirped QPM nonlinear crystals that involve nonlinear 
grating in the spatial coordinate along the direction of 
pump propagation with a nonuniform poling period [l7l - 
120( 1 . In this case, phase-matching is achieved along the 
crystal in such a way that the phase-matching conditions 
vary longitudinally by use of QPM which leads to broad- 
band biphoton. 

The theory of linearly chirped biphoton wave-packets 
is usually based on simple model allowing analytical de- 
scription of linearly chirped periodically-poled crystals 
or optical fibers (see, |17H20l |). In this model a phase- 
matching function of three- wave interaction exp[iAk(z)z] 
with a QPM grating with linearly varying spatial fre- 
quency Ak(z) = Ak — az has been considered. The 
interaction leads to biphoton spectral density in form 
of the erfi(x) error function. Such phenomenological 
approach qualitatively explains the basis properties of 
chirped quasi-phase-matching, in this number: extremely 
wide spectral bandwidths of biphoton wave packets and 
the sharp temporal separation between the signal and 
idler photons comprising a pair. However, this approach 
has been unable to describe some important details of 
chirped three-wave interactions. 

In this paper, we present more detailed description of 
chirped structures for generation of two-photon light via 
SPDC from a cw pump in another way: as the definite 
assembly of nonlinear segments. This approach provides 
more flexibility for designing of QPM gratings and al- 
lows us engineer the controllable phase relation between 
the various spectral components. On the whole it be- 
comes possible to control the frequencies and bandwidths 
of the signal and idler photons by varying also the num- 
ber of layers equally with the local poling period along 
the length of the crystal. In this approach the total am- 
plitude of two-photon generation is calculated as the su- 
perposition of partial amplitudes corresponding to each 
layer with local spatial chirp frequency. As frequencies of 
signal and idler photons are varied with different layers a 
much broader range of photon frequencies is formed due 
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to the superposition. Such detailed analysis of chirped 
QPM allows considering interference and transition ef- 
fects in SPDC stipulated by definite numbers of layers. 
In this approach the quasi- phase-matching function of 
SPDC is expressed through so called Gauss sums instead 
of the continuous error function that is appeared in the 
phenomenological approach. It should be mentioned that 
during recent years it has been seen an impressive num- 
ber of experiments implementing Gauss sums 



Sn(0 



M 

E 



w„ 



5 27ri(m+2SrK 



(1) 



-M 



in physical systems, particularly, for number factoriza- 
tion schemes. These systems range from nuclear mag- 
netic resonance (NMR) methods |2l| - |23j via cold atoms 
[IH and Bose-Einstein condensates (BECs) [25[ , tailored 
ultrashort laser pulses (2|| [27j to classical light in a multi- 
path Michelson interferometer [2^, [29[ (for an introduc- 
tion to this field, see [3(1 HU). In addition to these results 
here we present new quantum systems for implementa- 
tion of Gauss sums to factor numbers. 

We investigate two structures: chirped layered pho- 
tonic crystals and aperiodically poled layered crystals. 
In the phenomenological approach in which the phase- 
matching function is described by the continuous error 
function both structures give the same results distinct- 
ing only on the physical parameters. However, detailed 
description of these chirped structures shows important 
differences between them. 

The paper is organized as follows: In Sec. II, the brief 
description of SPDC in multilayered media is presented. 
Sec. Ill is devoted to chirp in layered photonic-like crys- 
tals. The production of broadband biphotons in aperiod- 
ically poled layered crystals is described in Sec. IV. We 
summarize our results in Sec. V. 



II. SPDC IN MULTILAYERED MEDIA: BRIEF 
DESCRIPTION 

In this section we briefly describe generation of two- 
photon light via SPDC in one-dimensional x*- 2 - 1 media 
consisting of layers with different coefficients of nonlin- 
earity and refractive indices [l(ifl . [l3| . fl4| . We consider 
collinear, type-II QPM configurations for generation of 
photons at the same frequency. Since the generated pho- 
tons are collinear, their directions cannot be used to dis- 
tinguish the two photons. Thus, we assume that photons 
in pairs have different polarizations, but omit the po- 
larization indexes below. Three-wave interaction Hamil- 
tonian is expressed as the sum of interactions in each 
layer in terms of the electric fields for the m-th layer, 
E 0m (z,t), EJ m {z,t),Ef m [z,t), j=(s,i), 



H(t) = 



Here x^ 2 H z ) 1S the second-order susceptibility, E 0m (z,t) 
represents classical laser field at the frequency u>q and 



dzx {2) (z)Eo m (z,t)E im( -) E { sm > (z,t) + h.c. 

(2) 



E S m(z,t) represent the positive-frequency 
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parts of the fields of the subharmonics centered at the 
frequencies u>i = ^ + fi, ui s = ^ — Q. 



The two-photon state can be written as 



duj s duji$(uj s , uJi)a + (iv s )a + (w, ) 1 0} , (3) 



where a + (uj s ) and a + (wj) are the creation photon opera- 
tors for modes with frequencies oj s and Wi, wo = w s 
0) is a vacuum state of the signal and the idler fields 
and $ (uj s ,LUi) is the spectral amplitude of two photon 
radiation. Two-photon amplitude for multilayered struc- 
tures has been presented in many papers [Tq|, [l3j], (see, 
also It has been calculated as the sum of defi- 

nite partial amplitudes in the general form given by the 
product of the pump envelope function Eq(ujo) and the 
phase matching function F(Ak) 
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iPm = l n Ak n , ¥>1 = 



(5) 



Here, Xm is the second-order susceptibility in the ru- 
th layer, l m = z m +\ — z m is the length and Ak m — 
fcg 1 — fc™ — k™ is the phase mismatch vector for the m-th 



layer, kf (z,uij) 



t (ui) is the correspond- 



ing refractive index of the medium at the given frequency 
that describes the effects of dispersion on the properties 
of photon pairs. In this case the probability of twin- 
photon generation is calculated as 



\®(Us,Wi)\ 
N 
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in the form which clearly demonstrates that the interfer- 
ing probability amplitudes A m ~ EoF m of two-photon 
generation in each layer contribute to the total probabil- 
ity. Further, for convenience we expand the wave vectors 
in Taylor series around the exact quasi-phase-matching 
frequency and take into consideration only the zero- and 
first-order terms: k s = k s (if)-Qk' s , h = h (^)+Qk' i , 
where k[ s are the first derivatives of the dispersion law 
evaluated at related to the group velocity. This ap- 
proximation means that the spectrum of biphotons is not 
too broad compared to the difference of group velocities 
of the subharmonics. 
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III. CHIRP IN PHOTONIC-LIKE CRYSTALS 

Recently, non linear photonic crystals have been used 
in the context of the process of spontaneous paramet- 
ric down-conversion. In particular, the semiconductor- 
based nonlinear one-dimensional photonic crystals have 

been explored for phase-matching; in generation of 

polarization-entangled photon pairs [32| , [33j . It has been 
demonstrated that one-dimensional photonic crystal are 
very effective for generation of photon pair due to field 
localization in such structures [HI, (3(|. The pro- 
duction of photon pairs with engineered spectral entan- 
glement properties in one-dimensional nonlinear photonic 
crystals has been shown in [37] . 

In this section, we study one-dimensional nonlinear 
photonic crystals for chirp configuration. In the scheme 
proposed generation of photon pairs by SPDC in \^ lay- 
ered structure characterized by a chirp in its linear opti- 
cal properties. Considering the general equation ([3]), we 
assume that layers have the equal lengths l m — I as well 
as the equal susceptibilities Xm — XOi while refractive 
indexes of the pump n(u)Q, z) — tiq(z) and the generated 
waves n(u>j,z) = rij(z) (J — s,i) have linear dependence 
on the coordinate and are given by 



n (z) -- 
rij(z) 



no, < z < I; 

uq — mlfto ml < z < (m + 1)/, 

rij, < z < I; 

rij — mlfij ml < z < (m + 1)1, 



(7) 



refractive indexes of the pump and subharmonic waves, 
correspondingly. In this case, the phase-mismatch vector 
for the m-th layer has a linear chirp of the following form 
Afc m — Ak — a(m — \)l, where Afc = fco(wo) — k s (uj s ) — 
ki{u>i) is the phase mismatch vector at the first layer, and 
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is the spatial chirp parameter. 



Taking into account ([7]) we present the general expres- 
sion ((5]) as 



N 



F(Ak) = l Xa e^ Fm, 



where sinc(x) 



sin(x) 



x sine 



( (Ak - a(m - 1)1)1 
V 2 



,(8) 



where /?o and j3j (j = s,i) are chirp parameters for the 



It is easy to see that in the frame of this presentation, 
the phase matching function is given by the Gauss sum, 
in spite of the error functions that are appeared in the 
phenomenological approach [1, [II]. Below, we calculate 
these sums considering the phase-mismatch function as 
Afc = Afc + MD, where Afc = k - k s (*f) - h(^) and 
D = k' s — k[ is temporal walkoff between signal and idler 
modes, assuming that Afco satisfies the QPM condition. 
In order to illustrate the broadening of the spectrum and 
the other features of obtained phase-matching function 
we have calculated the squared amplitude of the phase- 
matching function |F(f2)| 2 on the base of the formula ([6]) 
and the wave vector dispersion expansion as 



\F(n)\=i 2 Xo 
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f (Dn + Afc - a(m - 1)1)1 
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We present a detailed analysis of normalized bipho- 

I F 1 2 

ton spectral density |/| 2 = jt^i in dependence on the 

wavelength of signal (idler) mode. In this way, we se- 
lect in the Fig. ([T]) the results for typical values of 
pump wave length Ao = 0.458/im, the chirp parameter 
a = 1200cTO~ 2 , B = cD = 0.3 and crystal with the total 
length L = 0.8cm taken from the experimental results on 
LiTa0 3 [HI]. Note, that our approach allows to consider 
the general case of arbitrary number of layers on the base 
of the formula (9). To develop comparatively analysis we 
represent the squared amplitude of the phase matching 
function for the crystals with the fixed total length L, but 
with the various number of domains, particularly, for the 



cases of small and large numbers of layers. 

The results are depicted in FigsQIa,b) for N = 5 
and N = 10 layers and the phase-matching conditions 
Afco = Sal and Afco = bal, correspondingly. As we see, 
in these cases the biphoton spectra consist from definite 
separated spectral peaks corresponding to the number of 
domains. The frequencies of the spectral lines are given 
by Q m = [a(m — 1)1 — Afc ]D _1 according to the for- 
mula (9) provided that the spectral lines are resolved i.e. 
their widths are less than the intervals between them in 
the range of the wave-length broadening. Thus, there are 
critical ranges of domains number when the spectral lines 
are not resolved and chirp leads to the broadening of the 
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FIG. 1: Biphoton spectral density in dependence on the wave- 
lengths of signal (idler) field for the crystals with the length 
L = 8000/xm, dispersion parameter B = 0.3, chirp parameter 
a = 1.2 x 10 -5 fim~ 2 and different numbers of layers: N — 5, 
Afco = Sal (a); N = 10, Afco = 5al (b); N — 20, Afco = 10a/ 
(c) and N = 80 ,Afc = 40a/ (d). The cases of the chirp- 
ing parameter a = 6 x 10~ 6 pim -2 : N = 5, Afco = 3a/ (e); 
TV = 80, Afco = 40a/, (f). The red dashed lines correspond to 
the wavelength 2Ao = 0.916/im. 



biphoton spectra. 

If domains' number is increasing, the chirp is dis- 
played as a broadening of the biphoton spectra as it is 

depicted in FigsQJcd) for N = 20 and N = 80 la y~ 
ers and the phase-matching conditions Afco = 10a? and 
Afco = 40a/, respectively. One can see that the spectrum 
of the biphoton field is quite broad (from 900 nm to 1300 
nm wavelength ) and has nearly rectangular shape al- 
ready for N — 80 (see, Fig. QJd)). For further increasing 
of the number of domains the interference picture become 
smooth and the results qualitatively coincidence with the 
analogous results obtained on the base of simple model 

mm 

It should be noted that the length of biphoton spectral 
shape does not depend from the number of layers for fixed 



total length of nonlinear media but is only determined by 
the chirping parameter. We illustrate this statement on 
FigH] (e,f) by consideration of the other chirping param- 
eter a = 0.000006/iTO~ 2 that is two times less than the 
previous one, for N=5 and N=80 layers, correspondingly. 
As we see, increasing of the chirping parameter leads to 
the increase of the intervals between the spectral peaks 
for the case of small numbers of layers. For the case of 
large layers' numbers this increasing leads to the increase 
of the range of the wave-length spectral broadening. 
IV. CHIRP IN APERIODIC ALLY POLED 
CRYSTALS 



Recently, a aperiodic poling have been used for various 
applications. It has been provided not only for compen- 
sation of a natural phase mismatch but has been allowed 
to tailor the properties of emitted photon pairs using non- 
linear domains with variable lengths (chirped periodical 
poling). Domains of different lengths in an ordered struc- 
ture allow an efficient nonlinear interaction in an ultra- 
wide spectral region extending typically over several hun- 
dreds of nm [3-il [H, |39| . It has been shown that photon 
pairs generated in such structures can posses quantum 
temporal correlations at the timescale of fs. Other appli- 
cations of ultra- wideband biphotons include nonclassical 
metrology (40| and l arg e bandwidth quantum informa- 



tion processing [4l|, [4: 



We consider aperiodically poled crystals as a multilay- 
ered structure consisting of N layers with variation of 
lengths: the length of each layer is larger from the previ- 
ous by a chirp parameter £. So the length of n-th layer 
is given by the following expression l n = Iq + (n — 1)£, 
where Iq is the length of the first layer. We consider 
the quadratic nonlinearity Xn having a constant \ n — 
(_1)(™- 1 )^ value within each n-th layer . This struc- 
ture can be divided into N/2 domains such that each of 
them will be consist of two layers with reversed crystal 
axes. So the poling period A which is the length of the 
domain races by 2C,l for any next domain, which means 
it is dependent on the coordinate. Thus, it is easy to 
realize that the relation between chirp parameters ( and 

a can be written as follows C = —■ 

According to the formulas ([5]) we obtain two-photon 
spectral amplitude F(Ak) for the case of aperiodically 
poled crystal in the following form 



F(Afc) = |^exp 



. AfcZ 



m , 
£(-l) m exp - 



iAk [ mlo + 
I 



(to - 1) 2 C 



. ^ Afc(Z + (TO-l)C) 
sml 



(10) 



through the discrete Gauss sum, where Afc is the phase ton generation as a function of the frequency £1 is pro- 
mismatch function. In this case the probability of bipho- 
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2E E(-D f 



, silJ , p fi + Afc ) a ° + C(m 2 +P ' 1)) ) sin ({Da + Afa) (fo + C( ^" 1)) ) cos (p(Z>0 + Afc ) (z + C(2m +/ ' 2) 



(11) 



Here A&o = fcj ( ^ ) + k s ( ^ ) — fco is the phase mismatch vector at the central idler and signal frequencies tOi 
where ujq is the laser pump frequency. In this formula we have considered also Ak = Ako + Di L 

It is easy to realize that the biphoton spectra for two 
systems under consideration arc essentially different in 
form. The main difference is that biphoton spectra for 
aperiodically poled structures do not consist from sepa- 
rated spectral peaks in the case of small layers in opposite 
to the case of photonic crystal. If the chirp parameter £ 
equals to zero the spectrum (11) is reduced to the spectral 
function of biphotons for pure periodically poled config- 
uration that includes segments of length Iq with positive 
Xo and negative — xo susceptibilities that alternate one 
to the other 
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\F(n)\ 2 = i 2 oX l- 



sin* 



' N(Dn+Ak -f-)l 



( (Dn+Ak --?-)l 



(DO, + Ak )l 



(12) 



In the approximation N >> 1 this probability reads as 
D (£>n+Afc ) s 



F\n\ 2 ~ l 2 N 2 sinc 2 



n {DQ+Ak ) y We illustrate the pe _ 

culiarities of the biphoton spectra in dependence on Q for 
aperiodically poled structure as well as for periodically 
poled structure on Figs. 2 for N = 50 and Afco = 0. As 
we see, in the limit C — an d N >> 1 there occur only 
one narrow pick for positive O centered at ^ . Increasing 
of the chirp parameter broads the frequency spectrum. 

Typical biphoton spectra for aperiodic, chirped poling 
are shown in Figs. [3] in dependence on the signal-field 
(idler-field) wavelength (A s > 2Ao) for various layer num- 
bers and chirping parameters. For comparison of these 
spectra with analogous results obtained for the case of 
photonic crystals as well as with the experimental re- 
sults (see, for example, jTH) we use the parameters that 
are suitable for the LiTaC>3 aperiodically poled crystal, 
with B — 0.3, the whole length L = 0.8cm and for the 
laser pump frequency Ao = 0.458/im. 

Biphoton spectral density for an aperiodically poled 
crystal with layers number N=50, chirp parameter £ = 1 
and lo — 109.5, which corresponds to a smaller spatial 
chirp parameter a — 240cm -2 is shown in the Fig. [3ja). 
The analogous result for the case of crystal with the same 
number of layers, but with chirp parameter £ = 2.82 and 
Iq = 88.09, which corresponds to the spatial chirp pa- 
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FIG. 2: Biphoton spectra in dependence on Q frequency for 
non-chirped QPM grating (( — 0,1 — (a)) and chirped 

QPM grating (£ = 2.82, l = 88.09fim (b)). The parameters 
are: L = 8000/oti, A = 0.458/im, B = 0.3, N = 50. 



rameter a = 1200cm -2 is shown in Fig. [SJb). From 
comparison of these two results we can see that increas- 
ing of chirp parameter on two times lead to wide broad- 
ening of the signal- and idler-field spectra. Moreover, 
for each number of layers there is an optimal value of 
chirp parameter for which the spectrum is the broad- 
est. In Fig. IHb,c,d), we show the results of forming 
the biphoton spectra as number of layer increases pro- 
vided that the chirping parameter is fixed. For this goal 
we consider the number of layers: N = 50, N = 100 
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FIG. 3: Biphoton spectral densities in dependence on wave- 
length for aperiodically poled crystals with the following pa- 
rameters: N = 50, lo = 109.5/im, £ = 1/im, Ako = r 2 - 
(a); N = 50, l = 88.09^m, £ = 2.82^m, Afc = (b); 
N = 100, l = 52.225/wi, C = 0.55rtm, Ak = ^ (c); 
N = 160, lo = 35.51/xra, C = 0.18/im, Afc = f- (d^Here 
In = lo + n£) and B — 0.3, Ao = 0.458iim and red dashed 
lines correspond to the wavelength 2Ao 

and N = 160 and the chirping parameters are chosen so 
that again a — 1200cm -2 . As we can see from these pic- 
tures the spectra in these three cases have approximately 
the same form and they qualitatively match with the re- 
sult for the biphoton spectra broadening presented in [l9[ 
but are different in details. The wider broadening occurs 
for the higher number of layers: for the number of lay- 
ers N — 160 the spectral width achieves 0.2zim. Thus, 
we conclude that for the case of aperiodic chirping the 
width of spectrum is strongly dependent on the number 
of layers despite the refractive index frequency chirp con- 
sidered in the previous section. We also observe that in 



the case of aperiodically poled crystals the spectral lines 
corresponding definite segments are not resolved as it is 
demonstrated in the case of photonic crystals. 



V. CONCLUSION 

In conclusion, we have investigated production of 
biphoton wave packet in two chirped QPM structures 
presented as the ensemble of second-order nonlinear lay- 
ers. In this approach the total amplitude of two-photon 
generation has been calculated as the superposition of 
partial amplitudes corresponding to each layer with lo- 
cal spatial chirp frequency. In this way, the spectrum of 
spontaneously generated chirped photons has been calcu- 
lated through Gauss sums that involve realistic parame- 
ters of three- wave interaction in each layers. The detailed 
calculations have been done for two QPM structures: 
(a) chirped nonlinear one-dimensional multi-layered pho- 
tonic crystals; (b) aperiodic poled one-dimensional multi- 
layered crystals. The advantage of the presented ap- 
proach is that we have found out clear and easy cor- 
relation between number of layers the crystal consists 
of and the formation of broadband interference shape of 
the biphoton spectral rate. The physical mechanism of 
chirping for both structures have been also explained for 
physically realised parameters. This approach allows us 
to consider the case of small number of layers. Partic- 
ularly, it has been demonstrated that biphoton spectra 
for chirped photonic crystals with small number of lay- 
ers (N=5,10) consist from well-resolved spectral lines. It 
has been demonstrated two mechanisms for forming the 
broadband spectra of signal (idler) waves in SPDC for 
both configurations as number of layers increases and for 
various chirping parameters. 
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